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ABSTRACT 


An  equation  has  been  derived  describing  a  complete 
battery  discharge  for  the  case  when  the  current  density 
distribution  is  uniform.  The  battery  potential  during  dis¬ 
charge  is  given  as  a  function  of  time,  current  density,  polar¬ 
ization,  internal  resistance,  and  other  factors.  This  equa¬ 
tion  will  be  used  as  a  base  for  deriving  the  more  general 
case  where  the  current  density  over  the  face  of  the  elec¬ 
trodes  is  uneven.  This  equation  has  a  number  of  practical 
applications.  It  can  be  used  to  describe  battery  charges 
and  discharges,  capacities,  power  evolution,  and  predict 
capacities  on  the  basis  of  limited  data.  The  equation  can 
be  applied  to  the  charging  of  batteries  by  changing  the  signs 
of  some  of  the  terms  in  the  equation. 
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THEORETICAL  DESIGN  OF  PRIMARY  AND  SECONDARY  CELLS 
PART  III  -  BATTERY  DISCHARGE  EQUATION 


INTRODUCTION 

An  ultimate  goal  in  the  Laboratories  battery  program  is  to  determine  procedures 
for  designing  batteries  having  optimum  properties  such  as  minimum  weight  or  volume. 

The  derivation  of  a  battery  discharge  equation,,  which  gives  the  cell  potential  during  dis¬ 
charge  as  a  function  of  discharge  time,  current  density,  and  certain  other  factors,  has 
been  included  as  part  of  the  series  of  studies  to  reach  this  goal.  This  equation  can  be 
used  to  describe  battery  charging  and  discharging,  capacities,  to  predict  capacities,  and 
as  a  basis  for  deriving  the  equations  that  describe  the  discharge  of  a  battery  through 
constant  resistance. 

If  the  potential  drop  in  the  grids  of  a  battery  during  discharge  is  relatively  high  com¬ 
pared  to  the  potential  drop  in  the  electrolyte,  the  current  density  will  be  unevenly  dis¬ 
tributed  and  will  be  much  higher  at  the  top  than  at  the  bottom  of  the  electrode.  The  curr 
rent  density  distribution  over  the  face  of  the  electrode  as  well  as  other  cell  characteristics 
has  been  determined  mathematically  for  linear,  logarithmic,  cubic,  and  other  types  of 
polarization.*  An  example  of  what  can  occur  has  been  plotted  (Fig.  1)  where  the  ratio  of 
the  actual  current  density  to  the  average  current  density  is  shown  as  a  function  of  the 
fractional  electrode  height  where  0  represents  the  bottom  of  the  electrode,  0.5  is  the  mid¬ 
point,  and  1.0  is  the  top  of  the  electrode.  The  dotted  line  represents  the  case  where  the 
current  density  is  uniform  over  the  face  of  the  electrode.  This  ideal  condition  is  ap¬ 
proached  closely  by  most  batteries.  In  certain  cases  an  appreciable  deviation  from  this 
ideal  distribution  may  occur.  This  condition  is  most  often  encountered  in  tall  batteries 
or  in  batteries  where  optimum  characteristics  are  being  sought.  In  the  curve  in  Fig.  1 
the  current  density  at  the  top  of  the  electrode  is  four  times  the  average  current  density 
and  at  the  bottom  of  the  electrode  is  about  1/7  of  the  average  current  density.  This  con¬ 
dition  would  ’  Id  true  only  at  the  beginning  of  the  discharge,  whenever  we  are  dealing 
with  porous  i  •  odes,  because  the  high  current  density  at  the  top  of  the  electrode  causes 
the  active  mater iai  to  be  used  up  at  a  faster  rate  than  average.  The  uneven  distribution 
of  unused  active  material  that  results  would  affect  the  current  density  distribution  and 
tend  to  make  it  more  even.  If  optimum  cell  properties  are  to  be  obtained,  a  mathematical 
analysis  is  needed, of  the  cell  characteristics  throughout  a  complete  discharge.  One  of 
the  bases  for  such  an  analysis  is  given  in  this  report,  in  which  the  potential  of  a  battery 
during  discharge  is  determined  for  the  ideal  case  w'here  the  current  density  distribution 
over  the  face  of  the  electrode  is  uniform. 


MATHEMATICAL  ANALYSIS 
Basic  Assumptions 

The  mathematical  analysis  given  here  is  based  on  the  assumption  that  the  following 
conditions  are  applicable: 


*C.  M.  Shepherd,  "Theoretical  Design  of  Primary  and  Secondary  Cells.  Part  I  -  Effect 
of  Polarization  and  Internal  Resistance  on  Current  Density  Distribution,"  NRL  Report 
5211,  December  29,  1958. 
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Fig.  1  -  Example  of  uneven  current 
density  distribution 


1.  The  anode  and/or  cathode  have  porous 
active  materials. 

2.  The  electrodes  and  electrolyte  are 
homogeneous. 

3.  The  electrodes  are  parallel. 

4.  The  current  density  is  uniformly  dis¬ 
tributed  over  the  face  of  the  electrode,  as  it  is, 
approximately,  in  most  batteries. 

5.  The  distance  between  the  electrodes 
and  the  cell  case  is  small  compared  to  the 
width  of  the  electrodes  and  the  amount  that  the 
electrolyte  extends  above  the  electrodes  is 
small  compared  to  the  height  of  the  electrode. 
Under  these  conditions,  end  effects  on  the  cur¬ 
rent  density  distribution  are  negligible. 

6.  The  temperature  is  constant. 

7.  The  internal  resistance  is  constant. 

8.  The  battery  is  discharged  at  a  constant 
current. 

9.  The  potential  drop  during  discharge  can 
be  divided  into  three  parts:  (a)  a  potential  drop 
occurring  only  in  the  initial  portion  of  the  dis¬ 
charge,  (b)  a  potential  drop  due  to  internal 
resistance  j  and  (c)  a  potential  drop  due  to 
polarization. 


10.  The  polarization  is  linear  over  the  range  of  current  densities  covered  in 
the  battery  discharge. 


Definition  of  Terms 

t  =  the  time  at  any  point  during  the  discharge. 

i  =  the  apparent  current  density  in  amperes  per  unit  area. 

M  =  the  amount  of  active  material  per  unit  area  that  is  available  for  discharge  pur¬ 
poses  on  the  controlling  electrode  at  the  beginning  of  the  discharge’  It  is  not 
necessarily  the  total  amount  of  active  material  per  unit  area. 

Q  =  the  available  quantity  of  electricity  per  unit  area  and  can  be  calculated  by 
Faraday's  law  to  be  equj.1  to  kM  where  k  is  a  constant. 

Mj  =  the  amount  of  available  active  material  per  unit  area,  used  up  at  time  t .  Con¬ 
sequently  kMj  =  it. 

tj  =  the  capacity  of  the  battery  in  time  when  discharged  at  current  density  i. 

K  =  the  slope  of  the  polarization  curve. 
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L  =  the  internal  resistance  per  unit  area  of  the  cell. 

Wj  =  total  energy  evolved  at  time  t  expressed  in  watt  hours  per  unit  area. 

E  =  the  cell  potential  during  discharge. 

E,  =  a  base  potential  that  is  constant  in  value. 

E^  =  E,  -  Ki  -  Li  . 

Ef  =  the  voltage  at  the  end  of  the  discharge. 

kj  =  E,  -  Ef ,  a  constant. 

A,B  =  empirical  constants. 

=  the  active-material  current  density. 

The  active-material  current  density  at  any  time  t  during  a  discharge  is 
defined  by  the  two  following  conditions. 

1.  At  the  beginning  of  the  discharge  equals  the  apparent  current 
density  i. 

2.  At  any  time  t  during  the  discharge,  i,,^  is  directly  proportional  to 
the  amperes  per  gram  of  unused  active  material. 

Thus,  if  2/3  of  the  available  active  material  has  been  consumed  at  some 
point  during  the  discharge,  1/3  of  the  available  active  material  is  left  and 
is  equal  to  31  at  that  point. 

The  above  two  conditions  define  the  active -material  current  density  by 
the  equation 


(  “  'll 

(1) 

■^am  “ 

/  kM  \  . 

kM  -  kHij  ^ 

(2) 

am  “ 

(3) 

a 

3 

II 

(1  -  it/Q)'^i  . 

(4) 

Derivation  of  Equation 

If  all  factors  other  than  polarization  are  omitted  temporarily,  then  assumption  10  ' 
states  that  E,  the  cell  potential  during  discharge,  is  a  linear  function  of  the  active  mate¬ 
rial  current  density.  This  relationship  is  shown  in  Fig.  2  as  a  solid  line  over  the  range 
of  current  densities  that  give  a  linear'  polarization.  If  this  function  is  extrapolated  to 
^■in  =  0,  t'  R  intercept  is  E, .  If  the  slope  of  the  polarization  straight  line  is  K, 


am  ' 


E  =  E.  -  Ki 


(5) 
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Substituting  Eq,  (3)  in  Eq,  (5), 

The  term  Li  is  subtracted  to  take  care  of  the 
potential  drop  due  to  internal  resistance  and  Eq.  (6) 
becomes 

Fig.  2  -  Effect  of  linear 

polarization  on  potential  jg  evaluated  mathematically  a  set  of 

curves  is  obtained,  one  of  which  is  shown  in  Fig.  3  as  a 
solid  line  which  depicts  the  potential  E  as  a  function  of 
the  quantity  of  electricity  it .  However,  there  is  a  drop  in  potential  at  the  beginning  of  a 
battery  discharge  that  is  not  included  in.  Eq.  (7),  This  initial  drop  is  more  or  less  rapid 
and  is  shown  in  Fig.  3  as  a  dotted  line.  Consequently,  another  term  must  be  added  to 
correct  for  the  difference  in  potential  between  the  solid  line  potential,  calculated  from 
Eq.  (7),  and  the  dotted  line  that  represents  the  actual  discharge  potential.  It  has  been 
found  experimentally  that  the  expression  Ae'®  it/o^  where  A  and  B  are  empirical  con¬ 
stants,  gives  ah  evaluation  of  the  initial  potential  drop  that  is  adequate  in  most  situations. 
When  this  terni  is  added  to  Eq.  (7)  the  final  equation 

E  =  E,  -  K  i  -  Li  +  Ae^°  (8) 

is  obtained. 


Fig.  3  -  Battery  discharge  curve 


If  the  basic  assumptions  were  true  and  numerical  values  of  E^,  K,  Q,  and  L  were 
known,  Eq.  (8)  could  be  used  to  calculate  the  cell  potential  E  as  a  function  of  discharge 
time  t  at  various  current  densities  i  for  those  portions  of  the  curves  that  lie  to  the  right 
of  the  point  a  in  Fig.  3.  The  portion  of  the  curve  that' lies  to  the  left  of  point  a  is  gen¬ 
erally  of  less  importance.  However,  the  location  of  point  a  can  only  be  determined  from 
actual  discharge  data.  Moreover,  the  work  involved  in  determining  values  of  E,,  K,  Q, 
and  L  would  be  large.  Consequently,  it  would  be  easier  and  more  accurate  in  most  cases 
to  determine  discharge  curves  experimentally,  than  it  would  be  to  attempt  to  calculate 
them  by  use  of  Eq.  (8). 
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However,  Eq.  (8)  can  be  used  to  great  advantage  by  fitting  it  to  actual  discharge  data 
and  determining  numerically  the  empirical  values  of  E^,  K,  Q,  L,  A,  and  B.  The  result¬ 
ant  equation  describes  a  set  of  experimentally  determined  discharge  curves  made  at 
various  current  densities  by  giving  the  potential  E  as  a  function  of  the  time  t,  the  cur¬ 
rent  density  i,  and  the  numerically  calculated  values  of  E,,  K,  Q,  L,  A,  and  B.  Such  an 
equation  gives  an  accurate  description  of  the  battery  discharge  and  can  be  used  to  de¬ 
scribe  energy  evolution,  cell  capacity,  and  also  to  predict  cell  capacities. 


FITTING  THE  EQUATION  NUMERICALLY 
TO  EXPERIMENTAL  DATA 


There  are  a  good  many  approaches  that  can  be  used  to  determine  the  numerical 
values  of  Ej ,  K ,  Q ,  L ,  A ,  and  b  in  Eq.  (8)  fi’om  experimental  discharge  data.  Most  of 
^{hem  are  involved  and  time  consuming.  A  mOt.  •  xl  is  desired  that  is  simple,  easily  ap- 
/ plied,  and  gives  results  that  are  reasonably  ucc  ’.  to.  The  following  approach  was  the 
/  most  satisfactory  of  those  tested. 

The  first  case  to  be  studied  is  that  of.fitting  Eq.  (8)  to  a  single  experimentally  deter¬ 
mined  discharge  curve  shown  as  a  solid  line  in  Pig.  4.  Three  points  are  selected.  Point 
3  is  slightly  below  the  knee  of  the  curve.  Point  2  is  somewhat  to  the  left  of  the  knee. 
Point  1  is  farther  back  along  the  curve  but  not  so  far  back  that  it  lies  to  the  left  of  point 
a  and  is  caught  in  the  initial  part  of  the  curve  where  the  potential  is  dropping  rapidly. 


The  potential  at  point  1  is  indicated  by  substituting  Ej  and  itj  for  E  and  it,  respec¬ 
tively  in  Eq.  (8)  which  gives 

Ej  =  E,  -  K 

Similarly  for  points  2  and  3 

^2  =  E3  -  K 

E3  =  E,  -  K 

In  nearly  all  cases  the  numerical  value  of  Ae'®  decreases  very  rapidly  as  it  is 
increased  in  value,  and  as  a  consequence  its  value  at  points  1,  2,  and  3  becomes!  ex¬ 
tremely  small  and  is  considered  equal  to  zero.  ( 


Subtracting  Eq.  (10)  from  Eq,  (9)  gives 


Ej  -  Ej  =  Ki 


Qfitj  -  itj) 


(Q  -  itjlfO  -  iti) 
Similarly,  Eq.  (11)  is  subtracted  from  Eq.  (10) 

Qfitg  -  it;) 

(Q  -  itsXQ  -  itj) 
_  (it^  -  iti)(Q  -  itg) 

Ej  "  E3,  (itj  itjlfQ  —  itj) 


Ej  -  E3  =  Ki 


(12) 


(13) 


Dividing  Eq.  (12)  by  Eq.  (13)  gives 


El  “  Ej 


(14) 
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When  numerical  values  ol  E  and  it  are  taken  from  the  discharge  data  in  Fig.  4  and 
substituted  in  Eq,  (14),  a  linear  equation  in  Q  is  obtained  which  is  solved  to  give  a  nu¬ 
merical  value  of  Q.  By  substituting  this  value  of  Q  in  Eqs.  (12)  or  (13)  and  solving,  a 
numerical  value  for  K  is  obtained.  When  the  numerical  values  of  Q  and  K  are  substi¬ 
tuted  in  Eqs.  (9)  and  (10),  remembering  that  the  exponential  terms  are  essentially  zero, 
a  numerical  value  for  -  Li  is  obtained.  When  these  numerical  values  of  E,  -  Li,  K, 
and  Q  are  substituted  in  Eq.  (8),  omitting  the  final  exponential  term,  E  is  obtained  as  a 
specific  function  of  the  quantity  of  electricity  (it).  If  this  value  of  E  were  plotted  in  Fig. 
4,  it  would  follow  the  dotted  line  to  the  left  of  a,  and  to  the  right  of  a,  would  approximate 
the  solid  line  which  represents  the  battery  discharge.  The  difference  in  potential  between 
the  dotted  and  the  solid  line  in  most  cases  can  be  fitted  approximately  by  the  empirical 
term 


Equation  (8)  may  also  be  fitted  to  discharges  made  at  a  number  of  current  densities. 
A  moderately  high  current  density  curve  and  a  moderately  low  current  density  curve  are 
selected  from  the  available  discharge  curves  and  four  points  are  selected  on  these  curves 
as  shown  in  Fig.  5.  Points  3  and  4  are  taken  just  below  the  knees  of  the  curves  and  points 
1  and  2  are  taken  near  the  beginning  of  the  discharge  but  not  so  close  as  to  be  caught  in 
the  initial  potential  drop.  Under  these  conditions  the  value  of  A  e*®  ‘  is  negligible  and 
the  potential  at  point  2  is 

I 


E  -  K 


Q 


i,  -  Li, 


(15) 


where  Ej  and  i.tj  are  the  values  of  E  and  it . 


Fig.  4  -  Fitting  equation  to  a 
single  discharge  curve 


Fig.  5  -  Fitting  equation  to  discharge  curves 
made  at  various  current  densities 
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When  numerical  values  of  E,  i,  and  (it)  are  taken  from  the  discharge  data  in.  Fig.  5 
and  substituted  in  Eq,  (19),  a  quadratic  equation  in  Q  is  obtained  which  is  solved  to  give 
a  numerical  value  of  Q.  By  substituting  this  value  of  Q  in  Eqs.  (17)  or  (18)  and  solving, 
a  numerical  value  for  K  is  obtained.  When  the  numerical  values  of  Q  and  K  are  sub¬ 
stituted  in  Eq.  (15)  and  a  similar  equation  for  Ej,  two  simultaneous  equations  are  obtained 
which  may  be  solved  to  give  numerical  values  for  E^  and  L.  When  these  numerical  values 
of  E,,  L,  K,  and  Q  are  substituted  in  Eq.  (8),  E  is  obtained  as  a  specific  function  of  cur¬ 
rent  density  and  time  (it).  The  term  Ae'®  is  determined  by  difference. 

Experimental  discharge  data  from  Edison,  Ni-Cd,  silver-zinc-alkaline,  air 
depolarized-zinc-alkaline,  lead  acid,  and  dry  cell  batteries  have  been  plotted  as  dashed 
lines  in  Figs.  6-11. 


0  12  3  4  5  6 


It  (AMPERE-HOURS  PER  CELL) 

Fig.  6  -  Comparison  of  the  discharge  data  for  an  Edison  cell 
with  the  calculated  results  from  the  equation 

E=  1.320  -  0.00656  T  "p  *  ~  0.0799  i  +  0.228  e"  *  ’ 
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Fig.  7  -  Comparison  of  the  discharge  data  for  a 
sealed  Ni-Cd  cell  with  the  calculated  results 
from  the  equation 


il  (AMPERE -HOURS  PER  CELL) 

Fig.  8  -  Comparison  of  the  discharge  data  for 
a  silver  cell  with  the  calculated  results  from 
tb.e.  £  aquation  .  _  _ _ 

E  =  Either  1^1.831  +  0.02 

-  0.00514  (^-  o-:V27-Il)  ■  ^  °  ^ 

or 

1.S667  -  0.0004  (i  -  -o-.oWa  u)  ‘  ‘ 

whichever  is  larger  in  value. 


E  (VOLTS) 
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Fig.  9  -  Comparison  of  the  discharge  data  for  air 
depolarized  cells  with  the  calculated  results  from  the 
equation 
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Fig.  10  -  CompaAison  of  the  discharge  data  for  a  lead  acid  cell  with  the 
calculated  results  from  the  equation 
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Fig.  11  -  Coxnparison  of  the  discharge  data  for  No.  6  dry  CollfJ 
with  ^he  calculated  results  from  the  equation 

E  =  1.274  -  0.3343  (  ; - „  ;  )  i  +  0.  269  1  +  0,  36  e' ‘  ^ 

\ 1  -  0.039  it  I 


The  potential  E  in  volts  is  given  as  a  function  of  the  quantity  of  oloctl’i'’!*"' 
ampere  hours  per  cell.  The  current  density  is  given  in  amperes  per  coll,  i'ho  pollUU 
are  calculated  from  the  equations  that  were  fitted  numerically  to  each  set  of  dlscluirgc 
curves  by  the  method  just  described.  It  can  be  seen  in  each  case  that  there  Is  u  good 
agreement  between  the  experimental  data  and  the  fitted  equation  at  every  current  density, 
No  high  current  density  or  low  current  density  curves  were  omitted.  A  few  intermediate 
current-density  discharge  curves  were  deleted  in  Figs.  6  and  7  solely  to  maintain  visual 
clarity  in  the  graphs. 

Equation  (8)  is  being  fitted  to  discharge  data  for  a  large  number  of  various  types  of 
batteries  in  addition  to  those  shown  in  Figs.  6-11. 


EQUATION  DESCRIBING  BATTERY  CHARGE 

If  the  basic  assumptions  are  assumed  to  be  true  during  the  charge  of  a  battery,  then 
the  following  equation  would  hold  for  the  potential  of  a  battery  during  charge: 

E  =  +  K  i  +  Li  -  Ae'®  (20) 

Equation  (20)  can  be  obtained  from  Eq.  (8)  by  reversing  the  signs  of  the  last  three  terms. 

Equation  (20)  has  been  fitted  to  the  charge  of  a  Ni-Cd  cell  in  Fig.  12.  The  solid  lines 
represent  the  actual  data  and  the  points  are  calculated  from  the  equation. 


NAVAL  RESEARCH  LABORATORY 


11 


CALCULATED  POINTS 
O  I6.0AMP 
A  2.5AMP 
•  O.iSAMP 

THE  DASHED  LINES  ARE  EXPERIMENTAL 
CHARGE  DATA 


— -O-  _ _ _ _ a 

. _ _ _ _ — - _ — — 

. - 


_ o 


_ 0-' 


/ 

t'o 


A 


4  5  6  7  8  9 

it  (AMPERE-HCHJRS  PER  CELL) 


Fig.  12  -  Comparison  of  the  charge  data  of  a  Ni-Cd  cell 
with  the  calculated  results  from  the  equation 
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energy  evolved  during  discharge 

The  total  energy  that  has  been  evolved  from  the  beginning  of  the  discharge  up  to 
time  t  is  defined  by  the  equation 


t 


Substituting  Eq.  (8) 


i2  -  Li^  +  Aie-®  '*/® 

*  \Q  -  It; 

Wj  =  [e,  it+KQiln  (Q  -  it)  -  Li^t  -  AQ/Be-® 

Wt  =  E,it  -KQi  In  (1  -  it/Q)'^  -  Li^  t  +  AQ/B  (l  -  , 


(21) 

(22) 


(23) 


If  the  constants  from  the  equation  describing  the  discharge  of  the  Ni-Cd  cell  in  Fig. 

7  are  substituted  in  Eq.  (23),  the  total  energy  evolved  at  time  t  is  found  to  be 

Wt  =  1.25  it  -  0.0238  i,  in  (1  -  1.05  it)' ^  -  0.006i^  t  +  0.0248  (l  -  e ' ^ .  (24) 

In  Fig.  13,  Wj  in  watt  hours  per  cell  has  been  plotted  as  a  function  of  t  in  hours  for 
various  values  of  the  cell  current  in  amperes.  The  solid  lines  are  obtained  from  calcula¬ 
tions  based  on  numerical  integration  of  the  actual  discharge  data.  The  points  are  calcu¬ 
lated  from  Eq,  (24).  The  agreement  between  the  calculated  and  the  measured  results  in 
such  cases  depends  entirely  upon  how  well  Eq.  (8)  c^  be  fitted  to  the  original  discharge 
data. 
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Fig.  13  -  Calculated  vs  measured  energy  evolved 
during  the  discharge  of  a  Ni-Cd  cell 


CELL  capacity 

The  discharge  of  a  battery  is  considered  to  be  terminated  when  the  potential  has 
dropped  to  a  particular  value  which  is  called  the  cutoff  voltage.  A  constant  cutoff  voltage 
for  the  various  current  density  discharge  curves  can  be  used  satisfactorily  whenever  the 
potential  difference  between  the  curves  is  small,  as  they  are  in  Fig.  8. 

When  the  potential  difference  between  the  discharge  curves  is  large,  as  it  is  in  Fig, 

6,  a  constant  cutoff  voltage  would  not  be  satisfactory  since  a  value  low  enough  to  take  care 
of  the  6,0-ampere  curve,  would  be  entirely  too  low  for  the  0.5-ampere  curve.  It  would 
be  desirable  in  such  cases  to  have  a  cutoff  voltage  that  decreased  as  the  current  density 
increased.  This  condition  is  achieved  in  the  following  manner.  If  the  equation  for  the 
6.0-ampere  curve  in  Fig.  6  is  plotted,  omitting  the  final  exponential  term,  a  solid  curve 
is  obtained  that  approximates  the  actual  discharge  curve  to  the  right  of  a  value  of  it 
which  is  equal  to  approximately  2.0.  When  it  is  less  than  2.0  a  dotted  line  is  obtained 
which  intersects  it  =  o  at  E^.  Values  of  E,.  are  obtained  for  each  current  density  curve 
in  a  similar  manner.  The  cutoff  voltage  on  each  discharge  is  taken  to  be  some  suitable, 
constant  value,  kj,  such  as  0,25  volt,  below  the  value  of  E^  for  that  curve.  This  opera¬ 
tion  gives  a  value  for  the  cutoff  voltage  that  decreases  as  the  current  density  increases. 

E,,  may  be  determined  mathematically  by  substituting  it  =  0  in  Eq.  (8)  and  omitting  the 
exponential  term 

Ev  =  E,  -  Ki  -  Li  .  (25) 

At  the  end  of  the  discharge,  the  exponential  term  in  Eq.  (8)  is  negligible  in  value  and 
the  final  or  cutoff  voltage  is 

''  "  (q^)  ‘  w 

where  t  ^  is  the  capacity  in  time  at  the  end  of  the  discharge  for  the  particular  current 
density  i. 

Subtracting  Eq.  (26)  from  (25)  gives 


■f- 


(27) 
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IS 


t  -  *^2^ 

'*  K,i  ^  +  IC  j  i 

''2Q 

~  Ki  +  kj 


(28) 

(29) 


Equations  (28)  and  (29)  give  the  capacity  in  time  t  j  at  any  current  density  as  a  func¬ 
tion  of  K  and  Q  which  are  determined  by  fitting  Eq.  (8)  to  discharge  data. 

Capacities  calculated  from  Eq..  (28)  or  Eq.  (29),  have  been  plotted  in  Fig.  14  where 
the  time  in  hours  is  plotted  against  the  current  density  i.  The  Ag'  curve  was  taken 
from  the  first  step  in  the  discharge  of  the  cell  in  Fig.  6  and  the  Ag  curve  was  taken  from, 
the  second  step.  All  of  the  data  in  Fig.  14  was  taken  from  Figs.  6-11.  The.  axes  had  to 
be  shifted  in  most  cases  to  bring  all  the  curves  into  one  chart,  and  the  figures  in  the  x 
and  y  column  in  Fig.  14  are  the  quantities  that  must  be  multiplied  times  t,.  and  i,  re¬ 
spectively  to  give  the  true  values  of  the  capacities  and  the  current  densitir'-. 


Fig.  14  -  Comparison  of  calculated 
and  measured  capacities 


PREDICTION  OF  CELL  CAPACITIES 

Equation  (29)  is  a  function  of  the  two  quantities  K  and  Q.  It  was  shown  earlier  that 
numerical  values  of  K  and  Q  can  be  calculated  from  the  data  in  a  single  discharge  curve 
made  at  a  constant  current  density.  Numerical  values  of  E  and  (it)  for  three  points  on 
such  a  curve  are  substituted  in  Eqs.  (12)  and  (14)  and  solved  to  give  numerigal  values  of 
K  and  Q.  When  these  values  are  substituted  in  Eq.  (29),  it  can  then  be  solved  to  give  the 
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capacity  t ^  in  time  at  any  current  density  i,  Thus,  the  data  from  a  single  discharge 
curve  at  one  current  density  can  be  used  to  predict  capacities  at  other  current  densities. 

The  results  of  such  predictions  are  shown  in  Fig.  1 5  which  has  been  plotted  in  the 
same  manner  as  Fig.  14.  A  good  agreement  is  obtained  in  all  cases  between  the  pre¬ 
dicted  and  the  measured  capacities. 

A  single  current  density  curve  was  chosen  for  each  battery,  three  points  were 
selected  on  this  curve,  and  the  capacities  were  predicted  through  use  of  Eq.  (28)  or  (29). 

The  variation  in  the  predicted  capacities  was  found  to  be  small  when  the  positions  of 
the  three  selected  points  were  varied  or  when  a  new  discharge  curve  with  a  different 
current  density  was  used.  To  avoid  bias  in  the  presentation,  the  predicted  capacities 
shown  in  Fig.  15  are  those  obtained  from  the  first  calculation  that  was  made.  However, 
there  are  a  number  of  definite  limitations  to  the  use  of  this  method  for  predicting  capaci¬ 
ties.  It  is  essentially  an  e.xtrapoiation  technique  and  is  subject  to  the  errors  and  limita¬ 
tions  found  in  extrapolation.  Assume  that  all  the  true  capacities  were  known  and  that 
they  were  plotted  as  a  curve  on  log  paper  in  the  manner  shown  in  Fig.  15.  Very  serious 
errors  would  be  introduced  if  the  current  density  that  was  chosen  as  a  base  for  predict¬ 
ing  capacities  happened  to  lie  on  a  point  on  this  curve  where  the  slope  was  steeper  than 
approximately  -1.4.  In  the  case  of  the  dry  cell,  the  lowest  current  density  discharge 
curve  was  the  only  one  that  could  be  used  safely  as  a  base.  Deviations  may  also  occur  in 
the  predicted  capacities  at  current  densities  that  are  more  than  five  or  six  tirnes  that 
being  used  as  a  base  for  the  predictions.  The  five  highest  current  densities  on  the  dry 
cell  curve  in  Fig.  15  fall  outside  this  limit  and  it  can  be  seen  that  the  difference  between 
the  predicted  capacities  and  the  measured  capacities  for  each  of  these  current  densities 
is  well  above  the  average  for  the  remainder  of  the  data. 


i  (AMPS  PER  CELL) 


Fig,  15  -  Comparison  of  measured  cell  capacities 
with  capacities  predicted  from  a  single  discharge 
curve 
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A  more  thorough  mathematical  and  statistical,  analysis,  covering  a  wide  range  of 
batteries  and  discharge  conditions,  is  being  contemplated  in  an  attempt  to  determine  as 
accurately  as  possible  the  limitations  and  the  procedures  involved  in  predicting  cell 
capacities  by  this  technique. 


CONCLUSIONS 

The  potential  E  of  a  battery  during  discharge  at  constant  current  density  is  shown 
to  be 


E 


Li  +  A  e 


-B  it/Q 


where  E,  is  a  constant,  K'  is  the  slope  of  the  polarization  curve,  Q  is  the  available  quan¬ 
tity  of  electricity,  L  is  the  internal  resistance,  i  is  the  current  density,  t  is  time  and  A 
and  B  are  empirical  constants. 


This  equation  can  be  used  to  predict  battery  discharges  or  to  describe  them  by  fit¬ 
ting  it  to  available  discharge  data.  It  can  be  used  as  a  basis  for  deriving  other  equations 
which  describe  battery  charging,  characteristics  of  discharges  made  at  constant  resist¬ 
ance,  power  output,  and  capacities.  It  can  also  be  used  to  predict  capacities  over  a  wide 
range  of  current  densities  from  the  data  obtained  in  a  single  discharge  curve  made  at  a 
constant  current  density. 
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